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Magnetic and Structural Studies of the Quasi-Two-Dimensional Spin-Gap System 

(CuCl)LaNb207 

Makoto YOSHIDA^*, Nobuyuki Ogata\ Masashi Takigawa^^, Jun-ichi Yamaura\ Masaki Ichihara\ Taro 
KiTANO^, Hiroshi Kageyama^, Yoshitami Ajiro^, and Kazuyoshi Yoshimura^ 

^Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581 
^Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502 

We report magnetization, nuclear magnetic resonance (NMR), nuclear quadrupole resonance 
(NQR), and transmission electron microscopy (TEM) studies on the quasi-two-dimensional 
spin-gap system (CuCl)LaNb207, a possible candidate for the J1-J2 model on a square lattice. 
A sharp single NQR line is observed at the Cu and CI sites, indicating that both Cu and CI 
atoms occupy a unique site. However, the electric field gradient tensors at the Cu, CI, and 
La sites do not have axial symmetry. This is incompatible with the reported crystal structure. 
Thus the J1-J2 model has to be modified. We propose alternative two-dimensional dimer models 
based on the NMR, NQR, and TEM results. The value of the hyperfine coupling constant at 
the Cu sites indicates that the spin density is mainly on the d{'iz^ — r^) orbital {z \\ c). At 1.5 
K, Cu- and Nb-NMR signals disappear above the critical field Bc\ — 10.3 T determined from 
the onset of the magnetization, indicating a field-induced magnetic phase transition at Bd- 

KEYWORDS: NMR, NQR, TEM, (CuCI)LaNb207, spin-gap, quantum spin system 



1. Introduction 

Lovif-dimensional quantum spin systems with frustrat- 
ing interactions have recently been studied intensively 
both from experimental and theoretical aspects. The S — 
1/2 antiferromagnet on a square lattice with the nearest- 
neighbor and the next-nearest-neighbor interactions (the 
J1-J2 model) is a well known example. When the near- 
est neighbor interaction Ji is antiferromagnetic ( Ji > 0), 
the ground state has a Neel order with the wave vector 
(tt, tt) for J2/J1 < 0.38 while a collinear order with the 
wave vector (tt, 0) occurs for J2/ Ji > 0.52.^"^^ A spin liq- 
uid phase is proposed between these phases near J2/J1 
~ 0.5.^"'^^ A spin liquid phase is also proposed for fer- 
romagnetic Ji (Ji < 0) between the collinear phase for 
'h/Ji < —0.51 and the ferromagnetic phase for larger 
J2/J1.3) The series of compounds (CuAr)LaNb207 {X — 
CI, Br) is a candidate for the J1-J2 model with mixed 
ferro- and antiferromagnetic interactions. The struc- 
ture is reported to be tetragonal (space group PA/mmm) 
with the CuX square lattice^'*' as shown in Fig. 1. In 
fact, (CuCl)LaNb207 has a singlet ground state with a 
spin-gap^' while (CuBr)LaNb207 shows a collinear an- 
tiferromagnetic order. 

In (CuCl)LaNb207, the spin-gap is estimated to be 27 
K by fitting the susceptibility data to the isolated dimer 
model. The neutron inelastic scattering results also ex- 
hibit magnetic excitations with the gap of 2.3 meV, al- 
though the resolution is limited. On the other hand, the 
high-field magnetization measurements revealed a sud- 
den increase at the critical field Bd = 10.3 T.^' The 
Zeeman energy at Bd, gfJ-BBd/kB — 15 K, is appar- 
ently not sufficient to close the spin-gap at zero field. In 
addition, the magnetization saturates at Bc2 — 30.1 T.^^ 
The large interval between Bd and Bc2 indicates a large 
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Fig. 1. Schematic drawing of the crystal structure of 
(CuCl)LaNb207.^' 



width of the triplet dispersion. This is inconsistent with 
the neutron inelastic scattering data indicating that the 
triplet mode is almost dispersionless.^-* 

There are also open issues about the crystal structure. 
The neutron diffraction study indicates that the CI posi- 
tion deviates from the ideal (0, 0, 1/2) position with the 
C4-symmetry in the tetragonal space group PA/mmm.^^ 
If the structure deviates from the ideal square lattice, 
various ways to form a spin-gap become possible. ^-^ Un- 
ambiguous determination of the structure is crucial for 
good understanding of a frustrated quantum spin system. 

In this paper, we report magnetization, nuclear mag- 
netic resonance (NMR), nuclear quadrupole resonance 
(NQR), and transmission electron microscopy (TEM) 
studies on (CuCl)LaNb207. Our results indicate that 
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Fig. 2. X-ray diffraction patterns of the (a) aligned and (b) un- 
aligned powder samples. The inset illustrates the method for 
preparing the aligned sample. 



none of the Cu, CI, and La sites has the C4-symmetry, in- 
compatible with the reported tetragonal structure. Thus 
the J1-J2 model for the simple square lattice has to 
be modified. We propose two-dimensional dimer mod- 
els which arc consistent with our NMR, NQR, and TEM 
results. Cu- and Nb-NMR signals disappear above the 
critical field Bd — 10.3 T at low temperatures, indicat- 
ing a field- induced magnetic phase transition at Bd ■ 

2. Experiment 

The powder sample of (CuCl)LaNb207 was synthe- 
sized by the ion-exchange reaction of RbLaNb207 with 
CuCl2.^'^-' The magnetically aligned sample was pre- 
pared by rotating a cylindrical tube filled with the 
(CuCl)LaNb207 powder and epoxy (STYCAST1266) in 
a magnetic field B = 7 T. The magnetic field was applied 
perpendicular to the rotating axis as shown in the inset 
of Fig. 2. In this configuration, the magnetic hard axis 
tends to align along the rotating axis. The X-ray diffrac- 
tion pattern for the aligned sample was obtained with 
a MacScience M03X diffractometer using a monochro- 
matic Cu-Ka radiation as shown in Fig. 2 (a). Here 
the aligned axis was kept perpendicular to the scatter- 
ing wave vector. Compared with the diffraction pattern 
for the unaligned powder shown in Fig. 2 (b), the peaks 
corresponding to (OOn) planes are strongly enhanced in 
the aligned sample, where n is integer. This indicates 
that the c-axis is the magncitic hard axis. The rocking 
curve for these peaks with the half width of about 1.5° 
shows successful alignment. However, other peaks such as 
(102) or (110) are also observed in the aligned sample, 
indicating that a certain fraction of the sample remains 
unaligned. 

The magnetic susceptibility x was measured for the 
unaligned powder sample with a SQUID magnetome- 
ter (Quantum Design, MPMS). NQR measurements in 
zero magnetic field have been performed for the ^"^Cu 
and ^^Cl nuclei in the unaligned powder sample. Nuclear 
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Fig. 3. Temperature dependence of x measured for the unaligned 
powder sample at 1 T. The solid circles represent the raw data. 
The open circles represent x after subtracting tlic impurity con- 
tribution shown by the solid line. The inset shows the magnetic 
field (B) dependence of the magnetization at 2 K. The dashed 
line is the fit to the form M = NgSiXBBs{x) + AB. 



spin-lattice relaxation rates were measured by the inver- 
sion recovery method. NMR measurements in magnetic 
fields up to 11. 7T have been performed on the aligned 
sample with the magnetic field applied along the c-axis 
for 63.65cu, 35,37 93]^^^^ j^^^ 139lj^ nuclci. The TEM ex- 
periments were carried out at room temperature using a 
JEM2010F system with an operating voltage of 200 kV. 
The specimen was finely ground in methanol and then 
placed on a Cu microgrid mesh for TEM observation. 

3. Results cind Analysis 

We first analyze the temperature dependence of x in 
section 3.1. In section 3.2, we show the NQR results, 
which demonstrate that both Cu and CI atoms occupy a 
unique site and there is no apparent disorder. The elec- 
tric field gradient tensors at the Cu, CI, and La sites are 
investigated from the field swept NMR spectra in sec- 
tion 3.3. These results indicate deviation of the crystal 
structure from the reported one with a simple square lat- 
tice. In section 3.4, the temperature dependences of the 
magnetic hyperfine shifts at the Cu and CI sites are de- 
termined. Then the value of spin-gap is estimated. The 
hyperfine coupling constants are also determined at the 
Cu, CI, and Nb sites. In section 3.5, the high field be- 
havior above Bci is investigated at the Cu and Nb sites. 
Finally, we show the TEM result in order to provide ad- 
ditional information about the structure. 

3.1 Magnetic susceptibility measurements 

Figure 3 shows the temperature (T) dependence of x 
measured at the field of 1 T. The peak at 16 K and 
rapid decrease at lower temperatures are quite consistent 
with the previous report,^) indicating a singlet ground 
state. Below 4 K, however, x increases with decreas- 
ing temperature probably due to magnetic impurities. 
In order to investigate the impurity contribution, we 
also measured magnetization M up to 7 T. The inset 
of Fig. 3 shows the ^-dependence of M at 2 K. The 
data can be fitted by the sum of a term proportional 
to the Brillouin function Bs{x) and a B-linear compo- 
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Fig. 4. NQR spectra (a) at the ''^Cu sites (T 
at the sites (T = 2.9 K). 



15.1 



1.9 K) and (b) 



ncnt M ^ NgSBsix) + AB {x ^ gSfiBB/kBT). A good 
agreement was obtained for the values S = 1/2,A^ = 0.72 
%/Cu, and A^aMb^ = 3.0 x 10~^ emu/mol (TVa is the 
Avogadro constant) as shown by the dashed hnc in the 
inset of Fig. 3. The first term leads to the free spin con- 
tribution to X, NANgS^BBs{x)/B, due to isolated mag- 
netic impurities (the solid line in Fig. 3). By subtracting 
this term from the measured x, we obtained the corrected 
X shown by the open circles. The corrected x maintains 
a finite value (3x10""^ cmu/mol) as T goes to 0. It is 
also detected in the magnetization curve as the S-linear 
component. This is much larger than the typical value 
of the Van Vlcck susceptibility of Cu2+ (IQ-^ - 10"'' 
emu/mol). The NMR shifts at all sites do not show such 
large residual values at low temperatures as we discuss 
later. Hence the large residual susceptibility at low T is 
considered to be extrinsic, e.g. due to a secondary phase. 

3.2 NQR spectra and relaxation rates 

Figure 4 shows the NQR spectra at zero field for the 
^^Cu and ^^Cl sites, both having the nuclear spin I = 
3/2. A sharp single NQR line is observed for each nuclei. 
This shows that both Cu and CI atoms occupy a unique 
site and there is no substantial disorder for the Cu and 
CI positions. The NQR frequencies vnqr are determined 
to be 29.68 MHz for '^^Cu and 14.88 MHz for ^^G\. The 
values of vnqr and the spectrum shape show no tem- 
perature dependence from 10 K to 1.5 K, indicating that 
the structure is unchanged at low temperatures. 

Figure 5 shows the inverse temperature dependence of 
1/Ti at the ^^Cu and ^^■"'Cl sites at zero field. At both 
sites, 1/Ti is well fitted by an activation law 1/Ti oc 
exp(— E/fcsT). The activation energy is obtained as 
E/kB = 23 ± 2 K for Cu and E/Ub = 21 ± 1 K for 
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Fig. 5. Inverse temperature dependence of l/Ti at the ^"^Cu and 
^^Cl sites. The solid lines show the fit to an activation law 1/Ti oc 
exp(-E/fcBr). 




Fig. 6. Frequency- field diagram at 3 K at the ®^Cu sites in the 
applied magnetic fields B parallel to the c-axis. The dotted lines 
represent the resonance modes for the case of axially symmetric 
EFG around the c-axis. The solid lines are the fit to the results 

obtained from the Hamiltonian (2). The inset shows the Cu- 
NMR spectrum at 43 MHz. 



CI. Below 4 K, l/Ti deviates from the activation law 
and become independent of temperature below 3 K. The 
origin of the constant relaxation rate at low temperatures 
is likely to be impurities. 

3.3 Field swept NMR spectra and EFG tensors 

In the inset of Fig. 6, we show a typical field swept Cu- 
NMR spectrum for the aligned sample with the magnetic 
field B applied parallel to the c-axis. The spectrum has 
four sharp peaks assigned to the ^''Cu- and ^'"'Cu-NMR 
lines from the aligned part of the sample. The spectrum 
also shows a broad tail spreading to higher field. This 
tail probably comes from the unaligned portion of the 
sample. The resonance fields are determined precisely 
from the peak positions. The frequency-field (f-B) dia- 
gram can be constructed by plotting the resonance fields 
for various NMR frequencies between 26 and 55 MHz as 
shown in Fig. 6. 

The resonance frequency is generally determined from 
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Fig. 7. Frequency- field diagram at 3 K at the ^^Cl sites in the 
applied magnetic fields B parallel to the c-axis. The solid lines 
show the fit obtained from the Hamiltonian (2). The inset shows 
the Cl-NMR spectrum at 31 MHz. 
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Fig. 8. Field swept NMR spectrum at the '^^^La sites at 6 K and 
31 MHz in the magnetic field parallel to the c-axis. The inset 
shows the Sf/fo versus 1/i/g plot at 115 K. 



the following Hamiltonian including the magnetic hy- 
perfine interaction and the electric quadrupole interac- 
tion 



H 



h-fiB ■ I + 1 ■ K ■ B) 



+ 



eQ 



6/(2/ - 1) 



^ r 



a, 13 



(1) 



where h is the Planck's constant, Q is the nuclear 
quadrupole moment, 7 is the gyromagnetic ratio, K is 
the magnetic hyperfine shift tensor, 6a/3 is (a 7^ /?) or 1 
(a = and 14/3 = d^V/dad(i [a, (i = x, y, or z) is the 
component of the electric field gradient (EFG) tensor. 
If the Cu sites have the C4-symmetry around the c-axis 
in the tetragonal crystal structure, Vaf) is axially sym- 
metric, that is, Vafj = for a 7^ /3 and V^z = —^Vxx 
{Vxx = Vyy, z II c). If the magnetic field is applied 
parallel to the symmetric c-axis, three resonance modes 
appear at the frequencies and Ivnqr ± J^^rI, where 
UR = {1 + KcchB and unqr = SeQV,,/2I{2I - l)h, 
provided that all the components of K are much smaller 
than one. These results shown by the straight dotted 
lines in Fig. 6 are apparently in contradiction with the 
nonlinear behavior of the experimental data. This gives 
the direct evidence that the Cu sites do not have ax;ial 
symmetry around the c-axis. 

In the following, we assume that the c-axis is one of 
the principal axes of the EFG tensor at the Cu sites. 
This is the case, for example, when the Cu-Cl planes have 
the mirror symmetry. Although this assumption may not 
be strictly valid in the real structure, it is not possible 
to determine Va/s from the available data without any 
assumption. The nuclear Hamiltonian (1) is rewritten for 
the magnetic field parallel to the c-axis (2;-axis) as 



(2) 



H = - h^B{l + K)I, 

I' 

where K is the magnetic hyperfine shift along the c-axis 



and Pa = {3eO/2/(2/ - Although eq. (2) 

contains four parameters {vz, I'xi t^y, and K), the num- 
ber of independent parameters can he reduced to two 
{K and u^) from the relations i^x + i^y + i^z =0 and 
J^NQR = i/fc'l -l- (z^x — ^yY /^^ where the NQR frequency 
vnqr is determined to be 29.68 MHz in section 3.1. Four 
eigenvalues of this Hamiltonian can be calculated analyt- 
ically in the case of / = 3/2, resulting in six resonance 
modes. As shown in Fig. 6, the experimental data are 
well fitted by the solid lines obtained from eq. (2) with 



= 16.26 ± 0.02 MHz, 



13.37 ± 0.02 MHz, v„ 



-29.64 ± 0.02 MHz, and K ^ 0.44 ± 0.02 %. 

Figure 7 shows the f-B diagram at the '^^Cl sites ob- 
tained from the field swept NMR spectra of the aligned 
sample. A typical spectrum is shown in the inset of Fig. 7. 
The applied magnetic field B was parallel to the aligned 
direction (c-axis). The spectrum also shows the ""^^^La- 
NMR lines. The negative intensity is due to improper 
rf-pulse conditions, which were optimized for CI signal. 
The f-B diagram can be fitted by solving the Hamilto- 
nian (2) in the same way as has been done for the Cu 
sites. In Fig. 7, the solid lines show the results calculated 
for the parameter values = 14.17 ± 0.02 MHz, = 
-3.13 ± 0.02 MHz, Vy = -11.04 ± 0.02 MHz, and K = 
0.04 ± 0.02 %. 

Figure 8 shows the field swept ^"^^La-NMR spectrum 
at 6 K and 31 MHz in the magnetic field parallel to the c- 
axis. The seven peaks split by the quadrupole interaction 
for / =7/2 are clearly observed. Since the quadrupole 
interaction at the La sites is much smaller than that at 
the Cu and CI sites, the second order perturbation theory 
is sufficient to extract the values of i/q, and K. We assume 
the c-axis to be one of the principal axes for EFG. The 
resonance frequency between the states |m) and |m — 1) 
is obtained from the Hamiltonian (2) as 



{^x - Vyf 



12Z/0 



|m(m- 1) - — j 



(3) 
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Fig. 9. Temperature dependence of K at the ®^Cu sites for the 
magnetic field parallel to the c-axis. The dashed and dotted lines 

show the calculated results for the isolated dimer model and the 
two dimensional coupled dimer model, respectively. 
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Fig. 10. Temperature dependence of the field swept spectrum at 
the Nb sites at 75 MHz for the applied magnetic fields B parallel 
to the c-axis. 



where i^o = 7!?.^^^ If we define the frequency shift dvm = 
^'m^m-i — 1^0, the frequency shift of the central Hne 5ui/2 
obeys the relation, 



12 



(4) 



Therefore, Svij^/vq plotted against l/vl should yield a 
straight line with the slope (5/12) (i/^; — z/^)^ and the in- 
tercept K at l/i'o = 0. The inset of Fig. 8 shows the 
6v/vo versus I/vq plot at 115 K. We find that K at 
the La sites is almost zero (-8 x 10~^ ± 2 x 10~^ %) 
and \vx ~ Vy\ = 1.57 MHz. The value of is deter- 
mined from the frequency shift for the satellite lines as 
z/j; = \um^rn-i - i^-m+i«-m|/(2m - 1). The values of 
Vz, Vxi and Vy are determined to be 2.23 ± 0.02, —0.33 
± 0.02, and -1.90 ± 0.02 MHz, respectively. 

The values of Va determined by NMR are listed in Ta- 
bic I for the Cu, CI, and La sites. It should be noted that 
the values in Table I are determined on the assumption 
that the c-axis is one of the principal axes of EFG. The 
EFG tensor is not axially symmetric around the c-axis for 
all sites. That is, none of these sites has the C4-symmetry, 
which is incompatible with the reported tetragonal struc- 
ture. Thus the J1-J2 model for the simple square lattice 
has to be modified for proper description of this system. 
We discuss modified models in section 4.3. 

3.4 Hyperfine shift and coupling constant 

By repeating the above analysis for different temper- 
atures, we obtained the temperature dependence of K. 
Figure 9 shows the temperature dependence of K at the 
Cu sites in the magnetic field parallel to the c-axis. The 



Table I. Principal values of EFG tensors at the Cu, CI, and La 
sites. 



Cu 



CI 



La 



ly^ (MHz) 16.26 ± 0.02 
(MHz) 13.37 ± 0.02 
uy (MHz) -29.64 ± 0.02 



14.17 ± 0.02 2.23 ± 0.02 
-3.13 ± 0.02 -0.33 ± 0.02 
-11.04 ± 0.02 -1.90 ± 0.02 



shift consists of the contributions from the spin and Van 
Vleck susceptibility at the Cu sites, K = Kg+Kyy. Each 
term is related to the spin and the Van Vleck susceptibil- 
ity as Ks = {Ahf/NAHB)Xs- Kyv = iAorb/NAHB)xvv- 
The chemical shift is expected to be negligibly small com- 
pared with the spin and orbital shifts due to rf-electrons 
of Cu^+. As shown in Fig. 9, K decreases with decreas- 
ing T and is extrapolated to zero at T = 0. This means 
Kg and A^f are positive. Since Kyy is always positive, 
we can conclude that Kyy is nearly zero and Kg ap- 
proaches zero at T = 0. It should be noted that while 
X contains impurity contribution, K is free from impuri- 
ties, and thus provides direct evidence for a spin-gap. 

We estimate the spin-gap in two different models. The 
dashed line in Fig. 9 is the fit to the isolated dimer model 



K 



C 



1 



(5) 



from which wc obtained E' /k-g, = 21.8 K and Kyy = 
-0.02 ± 0.02 %. Note that the shift was determined in 
the range of magnetic field from 1 to 5 T. Considering the 
reduction of the gap by the average field gfJ.B{B)/k-B — 
4 K {{B) — 3T), the gap at zero field is estimated to 
be 26 K. This value is rather close to the gap (~27 K) 
estimated from the neutron inelastic scattering and the 
susceptibility measurements.*' However, the isolated 
dimer model ignores the dispersion of the triplet exci- 
tations. In the opposite limit, where the triplets have 
large dispersion Cq = E" + cq^, the spin susceptibility 
is determined by the excitation near the bottom of the 
dispersion Kg oc (l/T) / D{e)n{e)de oc e:^^{- E" / KbT) , 
where D{e) is the density of state, which is constant for 
e > E" in two dimension, and n(e) is the Bose factor. 
The dotted line in Fig. 9 shows the fit to this model with 
E" /Ub = 15 K. Therefore, the gap is estimated to be 19 
K at zero field. The gap values obtained in this study are 
compared with the reported results in section 4.1. 

Figure 10 shows the temperature dependence of the 
field swept NMR spectrum at the Nb sites. The nine 
peaks for / = 9/2 is observed at 1.5 K and the central 
line is indicated by the arrow. The central and the eight 
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Fig. 11. K-x plot for the Cu, CI, and Nb sites. 



satellite lines shift to lower fields with increasing tem- 
perature, indicating sizable hyperfine coupling at the Nt 
sites to spin magnetization of Cu ions. In addition, we 
observed another peak marked by the solid square ir 
Fig. 10. This peak does not shift with temperature, in- 
dicating no hyperfine coupling to the Cu spins. In fact, 
1/Ti for this peak is ten times smaller than that for the 
intrinsic signal at low temperatures. Therefore, we con- 
clude that this signal comes from extrinsic nonmagnetie 
phases. Imperfection of the ion-exchange reaction is a 
possible origin of such nonmagnetic impurely phases. 

Figure 11 shows the K-x plot for various sites obtaineej 
for the temperature range 1.5 ^ 20 K. Here, a Curie 
term due to magnetic impurities was removed from the 
X data as described in section 3.1. The hyperfine coupling 
constants Ahf for B \\ c are determined to be 13.8 ± 0.3 
6.9 ± 0.5, and 2.84 ± 0.08 T//zb for the Cu, CI, and Nb 
sites, respectively, (Table II). Although the NMR shifts 
at all sites approach near zero as T ^ 0, x maintains a 
finite value (3 x 10"'^ emu/mol) as mentioned in section 
3.1. Note that the impurity term from free spins has been 
already subtracted. The residual x a-t low temperatures, 
which corresponds to the iJ-linear magnetization in the 
inset of Fig. 3, does not couple to any nuclei observed by 
NMR. Hence we conclude that it has as extrinsic origin. 

3.5 High magnetic field region 

Figure 12 shows the field swept NMR spectrum at 
1.5 K for the Cu sites at different frequencies in the 
higher magnetic field region above Bci — 10.3 T. The 
very sharp peaks are due to the NMR-coil made from 
Cu. The ^'^'^^Cu-NMR lines from the aligned sample are 
indicated by the arrows. Their intensity decreases with 
increasing B and the peaks from the aligned sample dis- 
appear slightly below B^ =10.3 T. This is because the 
spin-spin relaxation time T2 is getting shorter when the 
field approaches Bd ■ There remains a broad signal near 



Table II. Hyperfine coupling constants for B \\ c. 



Fig. 12. Frequency dependence of the field swept NMR spectrum 
at the Cu sites at 1.5 K in the applied magnetic fields B parallel 
to the c-axis. The dotted lines show the critical field Bd ■ 
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Fig. 13. Frequency dependence of the field swept NMR spectrum 
at the Nb sites at 1.5 K in the applied magnetic fields B pajrallel 
to the c-axis. The inset shows the BcenteT-dependence of the 
integrate intensity. The dotted lines show the critical field Bd ■ 



Bel, observed in the spectrum at 130 MHz. Since this 
signal has much longer T2, it is probably due to some 
impurity phases. 

Figure 13 shows the frequency dependence of the field 
swept spectrum at the Nb sites. The quadruple split 
peaks for ^^Nb (/ = 9/2) are observed at 90 MHz. The 
intensity of the split peaks decreases with increasing B 
and the signal disappears slightly below Bci in a simi- 
lar way as has been observed for the Cu signals. Above 
Bel, the spectrum with one central peak and broadened 
quadrupole satellites is observed (115 MHz). This sig- 
nal probably has the same origin as the unshifted signal 
observed at 75 MHz at high temperatures (Fig. 10). In 
the inset of Fig. 13, the integrated intensity of the whole 
spectrum is plotted against the field of the central line 
(Bcenter)- The intensity normalized by the value at 7.2 T 
decreases with increasing Bcenter and vanishes at 9.5 T, 
slightly below Bd, except for the contribution from im- 
purity phases, indicating that a field-induced magnetic 
phase transition occurs around Bci- The field-induced 
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Fig. 14. (a) High-resolution TEM image of (CuCl)LaNb207 at 
room temperature along the [001]-zone axis and (b) its enlarged 
view, (c) Corresponding electron diffraction pattern, where in- 
dices are given based on the tetragonal unit cell of the PA^jmrnm 
structure.^) 



transition is also indicated by the recent specific heat 
measurements in magnetic fields. 

3.6 TEM result 

In order to obtain further information about the crys- 
tal structure, we performed TEM measurements at room 
temperature. Shown in Fig. 14 (a) and (b) are the [001] 
zone lattice images that demonstrate contrasts with a 
special period of about 4 A along (100), consistent with 
the lattice parameter. However, the corresponding elec- 
tron diffraction (Fig. 14 (c)) has revealed weak reflec- 
tions such as (1/2 0) or (0 1/2 0), in addition to strong 
fundamental reflections relevant to the reported crystal 
structure. The weak reflections are commensurate and in- 
dicate the doubling of the lattice period of both a- and b- 
axes, i.e., the size of the unit cell is given by 2a x 2b. Such 
a superstructure has not been detected even by high q- 
resolution synchrotron X-ray diffraction,^'^) presumably 
due to tiny modulation of the atoms and/or small domain 
size. Note that the superstructure reflection at (1/2 1/2 



0) was observed in (FeCl)LaNb20^ 
4. Discussion 



13) 



4.1 Spin- gap 

The gap value in (CuCl)LaNb2 07 is still controversial, 
while it has been investigated by several experimental 
methods. The most direct method is the inelastic neutron 
scattering, which shows the magnetic excitation with the 
gap energy of 2.3 meV (26.7 K).*) Although it is possible 
that the limited energy resolution masks fine structure of 
the neutron spectrum and the true gap may be somewhat 
smaller, the large discrepancy between this value and the 
gap derived from the critical magnetic field, g/XB-Bci/^B 
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Fig. 15. Orbital energy levels for Cu^"*" in a tetragonal crystal 
field. The ground state wave function is dix^ ~ V^) ™ or 
d(3z2 -r2) in (b). 



= 15 K, cannot be explained by such an effect. 

Prom the susceptibility measurements, the gap is esti- 
mated to be 27 K by using the isolated dimer model.'*) 
The magnetic hyperfine shift determined by NMR is free 
from impurity contribution and provides more reliable 
spin susceptibility. Using the isolated dimer model, the 
gap value is determined to be 26 K, consistent with 
the susceptibility results. We have also examined a two- 
dimensional coupled dimer model with a large width of 
the triplet dispersion and obtained the gap value of 19 K. 
Since these two models correspond to the opposite limits 
of zero and infinite width of the dispersion, the values 26 
K and 19 K set the upper and the lower bounds for the 
true gap. From the NQR relaxation rate 1/Ti at the Cu 
and CI sites, the gap value is determined to be 23 ± 2 and 
21 ± 1 K, respectively, consistent with the above range. 
This range is still significantly larger than gjjL^Bc\/k-Q (= 
15 K). 

To reconcile the discrepancy between the gap esti- 
mated at low magnetic fields and Bc\, Kageyama et al. 
proposed formation of two-triplet bound states with S 
= 2.^) If two-triplet bound states are stabilized due to 
frustration, the gap can close before B reaches A/gfXB 
and a nematic order may be realized above Bd-^^^ The 
Cu- and Nb-NMR signals were undetectable above Bd 
at 1.5 K due to divergence of I/T2 as mentioned before 
(Figs. 12 and 13). To investigate the magnetic structure 
above Bci, we need measurements at higher fields and 
lower temperatures. 

4.2 d-orbital 

The experimental determination of Ahf and Kyy at 
the Cu sites provides information about the Cu d-orbital. 
In this section, we compare the experimental results of 
Afif and Kyy with the calculated values. Figure 15 shows 
the orbital energy levels for Cu^+ in tetragonal crystal 
fields. The Cu'^^ ions have one d-hole with either d{x'^ — 
y^) or d(32^ — r^) wave function. The hyperfine coupling 
constants (for B \\ z) and A^ (for B ± z) are 
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expressed as 



Q 
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(r-3)(-« 
{r-')[-n 



An 



11 A 



8A 



2 

+ 7+ 7A 

4 6 A 



1^ 



45 A 



for d{2z^ -r^), 
(6) 

where /io is the permeabiUty of vacuum, k is the core 
polarization coefficient, A is the spin-orbit couphng pa- 
rameter, and Ai (i = 0, 1, 2) is the energy between the 
ground state and the excited sate shown in Fig. IS.""^^'"^^^ 
The values of k. A, and {r~^) for Cu'^+ ions are common 
to all insulating materials within 10-20 %, typically given 
as K = 0.28, A = -8.8 x 10"^ eV, and (r'^) = 4.26 x lO^i 
15,16) depends on materials but should be of the 
order of a few eV. We take A^ — 2.0 eV.^^'^^' By using 
these values, and A^^^ are estimated to be —24: A and 
2.71 T//XB for rf(.x2 - y"^) and 10.2 and -12.2 T//iB for 
d{3z^ —r"^), respectively. Compared with the experimen- 
tal values, Ahf (=13.8 T/^b), Ay for d(a;2 - y'^) and 

aY for diSz"^ — r^) can be excluded because they lead 

to large negative shifts. Note that the magnetic field is 
parallel to the c-axis in the experiment, and we consider 
the cases where the 2:-axis of the orbital is either paral- 
lel or perpendicular to the c-axis. Among the other two 
cases, AY for d{3z'^ — r^) is by far the favorable. 

The vanishing Van Vleck susceptibility Kyy — for 
B \\ c also supports the (i(3z^ — r^) orbital. For B \\ z, the 
Van Vleck susceptibility xvv (emu/mol) is calculated as 

\{^n\LM\' 



Xvv = A^A/Ue 



En — eo 



(7) 



where ipo is the grand state with the energy eo and 
^pn is the excited state with the energy e„. For the 
d{'iz^ — r^) orbital, xvv is zero. On the other hand, xvv 



for d{z 



is calculated to be A^aMb/Ai — 1-6 x 10 



cnm/mol, leading to Kyy = 2{r-^)xvv/NA = 0.23 %. 
This is clearly incompatible with the experimental ob- 
servation. Therefore, we conclude that the Cu spins are 
on the rf-orbital consisting mainly of d{3z^ — r^). This 
is natural because the Cu-0 bonding (1.84 A) is much 
shorter than the Cu-Cl bonding in the ab plane (3.14 or 
2.40 A).^) Such structure and the difference of the ionic 
valence (0^~ vs. Cl~) make the d{3z^ — r^) orbital fa- 
vorable for the rf-hole. It should be noted that broken 
axial symmetry revealed by our NMR data should allow 
slight mixing of other components such as d{x'^ ~ V^)- 
The large hyperfine field at the Nb sites can then be 
understood as due to the transferred hyperfine interac- 
tion through the d{2>z'^ — r^) orbital. We suspect that the 
interplane exchange interaction through the path Cu-0- 
Nb-O-Nb-O-Cu and the nearest-neighbor exchange in- 
teraction through the path Cu-O-O-O-Cu may not be 
completely negligible. 




\ — 
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Fig. 16. Possible structures for the CuCl plane. The open and 
solid circles represent the Cu and CI ions, respectively. The 
dashed lines in (a), (b), (e), and (f) indicate the unit cell for 
the superstructure. In the structure (e) and (f), the solid and 
gray spheres represent the CI atoms above and below the Cu 
plane, respectively. 



4-3 Possible model 

Precise structural information of the CuCl plane is cru- 
cial to understand the origin of the spin-gap in this com- 
pound. Our NQR results show that both Cu and CI have 
only one site and there is no apparent disorder. This 
means that periodic superstructure should be realized, 
even though the domain size of the superstructure may 
be very small. Although Caruntu et al. proposed disor- 
der in the CI position based on their neutron diffraction 
experiments,^^ their results are not necessarily incom- 
patible with the NQR results because possibility of a 
superstructure was not considered in the analysis of the 
neutron data. In addition to the NQR results, the EFG 
tensor does not have axial symmetry at the Cu, CI, and 
La sites and the TEM result shows that the unit cell is 
doubled along both the a- and the 6-axes from the origi- 
nal square lattice. Therefore, J1-J2 model for the simple 
square lattice has to be modified. Whangbo and Dai pro- 
posed models of ring clusters with even number of spins 
based on the neutron results to account for the spin-gap. 
They lead, however, to several inequivalent Cu and CI 
sites, inconsistent with our NQR results. 

The rather large spin-gap suggests that the real struc- 
ture has such a feature that allows natural formation of 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 9 



Cu dimcrs. Figures 16 (a), (b), (e), and (f) show examples 
of such structures, where we assume that Cu maintains 
the ideal position for the square lattice, (1/2, 1/2, 1/2). 
The structure in Fig. 16 (a) is simplest and leads to a 
two-dimensional dimer system as shown in Fig. 16 (c). 
Here, dimers are formed by the dominant exchange path 
Ji, and connected with each other by J(, J{', and J2. If 
Ji is antiferromagnetic, a spin-gap is naturally formed. 
This structure has only one Cu and CI site. The neutron 
diffraction study suggests that CI moves off the ideal (0, 
0, 1/2) position to one of the {x, 0, 1/2), {-x, 0, 1/2), (0, 
X, 1/2), and (0, —x, 1/2) positions in a random manner. 
The structure in (a) is an ordered version of the struc- 
ture proposed from the neutron experiments. However, 
the unit vectors a ± 6 of the structure (a) leads to the 
superstructure reflection at (1/2 1/2 0) but cannot ac- 
count for the observed TEM reflections such as (1/2 
0) and (0 1/2 0). Therefore, CI should take more gen- 
eral positions. Wo propose three models (b), (e), and (f), 
which are consistent with both the NQR and TEM re- 
sults. In the structure (b), CI moves off the ideal (0, 0, 
1/2) position to a position, (x, y, 1/2), {—x, y, 1/2), {x, 
—y, 1/2), or {—X, —y, 1/2), still maintaining only one 
site. The dominant exchange path Ji makes the dimers, 
which are connected by other exchange paths J[, Jf, J2, 
and J2 as shown in Fig. 16 (d). In the structures (c) and 
(f), CI moves off to {x, 0, 1/2+S), (x, 0, 1/2-6), {-x, 0, 
1/2+S), or (—X, 0, 1/2— S). The solid and gray spheres 
represent the CI atoms above and below the Cu plane, 
respectively. In both (e) and (f), the exchange paths are 
the same as the structure (a) that is shown in Fig. 16 
(c). Both models (c) and (d) generate frustration. 

A useful insight is obtained by comparing the magnetic 
ground state of (CuCl)LaNb207 and (CuBr)LaNb207, 
which shows a collinear antiferromagnetic order. We 
have observed a large internal field of 16 T at the '^^'^^Br 
sites in (CuBr)LaNb207 by zero field NMR. The de- 
tails will be published in a separate paper. If Ji is an- 
tiferromagnetic in (CuBr)LaNb207, the internal field at 
the Br sites should be canceled out for the structures 
(a), (e), and (f). Therefore, if the two materials have 
the same structure, (b) is the only possibility. However, 
the exchange interaction through CI or Br is very sensi- 
tive to small structural variation and often changes from 
ferromagnetic to antiferromagnetic.^^' Therefore, Ji for 
(CuBr)LaNb207 may be ferromagnetic, even if ,/i for 
(CuCl)LaNb207 is antiferromagnetic. If Ji is ferromag- 
netic in the structures (a), (e), and (f), a large internal 
field of the ordered moment survives. Therefore, both 
models (c) and (d) are possible. It seems that these ex- 
amples exhaust possible structures within the condition 
of the uniqueness of the Cu and CI sites and the doubled 
periodicity along the a- and 6- axes, although we are not 
able to prove rigorously. 



5. Summary 

NMR, NQR, and TEM results of (CuCl)LaNb207 have 
been reported. A sharp single NQR line observed at each 
Cu and CI site indicates that both Cu and CI atoms 
occupy a unique site. The electric field gradient tensor 
shows that there is no axial symmetry around the c- 
axis at the Cu, CI, and La sites. These results indicate 
that the J1-J2 model for the simple square lattice has to 
be modified. We propose a two-dimensional dimer model 
for this compound. Cu- and Nb-NMR signals disappear 
above the critical field Bd ~ 10 T at low temperatures, 
indicating a field-induced magnetic phase transition oc- 
curs at Bel- 
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